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Effects of medium polarization are studied for So pairing in neutron and nuclear matter. The 
screening potential is calculated in the RPA limit, suitably renormalized to cure the low density me- 
chanical instability of nuclear matter. The selfenergy corrections are consistently included resulting 
in a strong depletion of the Fermi surface. All medium effects are calculated based on the Brueckner 
theory. The 1 So gap is determined from the generalized gap equation. The selfenergy corrections 
always lead to a quenching of the gap, which is enhanced by the screening effect of the pairing po- 
tential in neutron matter, whereas it is almost completely compensated by the antiscreening effect 
in nuclear matter. 
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I. INTRODUCTION 

A satisfactory description of superfluidity in nuclear 
matter has not yet been achieved despite almost fifty 
years of research have elapsed since the first application 
of the BCS theory to nuclear systems Q. Somewhat 
at variance with the electron pairing in superconduc- 
tors the pairing in nuclear systems results from the in- 
terplay between the direct action of the bare nuclear 
force and the action induced by the medium polariza- 
tion. The attractive components of the bare nuclear in- 
teraction have led to the investigation of several pairing 
configurations, e.g. neutron-neutron or proton-proton 
pairing in the Sq channel in neutron stars _2j disre- 
garding possible repulsive effect exerted by screening of 
the force via the medium. A pairing suppression has in 
fact been found by most calculations of pairing in neu- 
tron matter (see Ref. and references therein). On 
the contrary, other pairing configurations have not yet 
been explored since the repulsive components of the di- 
rect nuclear interaction cannot support the formation of 
Cooper pairs. But there are strong indications that, in 
a nuclear rather than neutron matter environment, the 
medium polarization of the interaction can favor the for- 
mation of Cooper pairs similar to the lattice vibrations 
in ordinary superconductors. These indications come 
both from nuclear matter calculations and from finite 
nuclei. In nuclear matter the medium enhancement of 
neutron-neutron 1 Sq pairing is to be traced back to the 
proton particle-hole excitations 0, and in finite nuclei 
to the surface vibrations Q- 

Another distinctive feature of the nuclear environ- 
ment is the presence of strong short range correlations 
that induce two effects relevant for the pairing: one is 
the depletion of the Fermi surface, which is experimen- 
tally supported by measurements of electron scattering 
on 208 Pb 0, the other one is the strong mass renormal- 
ization caused by short-range particle-particle correla- 
tions pfl. 

The two effect conspire against the pairing formation: 



The depletion of the Fermi sea reduces the phase space 
available for particle-particle virtual transitions around 
the Fermi surface, the mass renormalization enhances 
the dispersive effect of the mean field ||- 

Therefore a complete microscopic treatment of the 
very subtle pairing problem requires vertex and selfen- 
ergy corrections to be treated and to be considered on 
the same footing. In a previous paper Q we made a 
study of these in-medium effects under several simplify- 
ing assumptions. First came the approximation to re- 
place the Born term of the pair interaction in the S=0, 
T=I channel by the Gogny force [jj. Though in that 
channel the Gogny force is not dissimilar to the action 
of the bare force ( see Ref. [lO|), it shows a little too 
much attraction for momenta characterizing saturation. 
The first improvement in the present work is then the 
use of a realistic two body force (VI8 see below) 
in the Born term. Secondly we will use as vertices in 
the induced interaction a force which is based on a more 
modern G-matrix calculation as this was the case 
for the Gogny force 9]. Thirdly we corrected an un- 
fortunate phase error which slipped into the evaluation 
of the induced force at least for the symmetric nuclear 
matter case in the S=0, T=l channel what gave raise 
to a too strong anti-screening effect. With these correc- 
tions and improvements we now get reasonable renor- 
malization effects of the pairing force and we calculate 
the corresponding gaps as a function of density in pure 
neutron matter as well as in symmetric nuclear matter. 

In detail the paper is organized as follows. In Sec. 
II the generalized gap equation is reviewed along with 
the approximations on the pairing potential and self- 
energy, which lead to the determination of the energy 
gap. In Sec. Ill the screening interaction is discussed 
in the RPA limit, and then the summation of bubble 
diagrams and the resummation of dressed bubble dia- 
grams both using the Landau parameters are derived. 
In Sec. IV the results are presented: first, for the sep- 
arate contributions of particle-hole (ph) scalar, vector 
and isovector excitations in neutron matter and nuclear 
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matter; second, the solution of the gap equation for Sq 
pairing with a discussion of the effects of the selfenergy 
corrections and medium polarization potential. Section 
V is devoted to the comparison with other calculations 
and to the conclusions. 



II. GENERALIZED GAP EQUATION 

The spectrum of a superfluid homogeneous Fermi sys- 
tem is derived from the generalized gap eauation |l3tll4l 



(1) 



where V is the sum of all irreducible NN interaction 
diagrams and F k (co) is the anomalous propagator. The 
class of diagrams selected for the present calculation is 
plotted in Fig. ^ 

In nuclear matter V can be approximated by the bare 
interaction Vo (diagram (a)), which is responsible of the 
formation of Cooper pairs, and the class of bubble in- 
sertions, which play the role of screening. In turn, the 
screening interaction Vi can be split into two parts: 
the one-bubble term (diagram (b)) containing only the 
mixed configuration with particle-particle (pp) plus ph 
excitations, the multi-bubble term also containing all 
insertions of pure ph interaction vertices (diagram (c)). 
The splitting is a convenient way to point out that 
mixed vertices and pure ph vertices have to be treated 
on different footing, as discussed in Ref. 4]. The first 
bubble diagram can also be seen as the lowest order cor- 
rection to the Born term. As discussed afterwards, the 
vertex insertions in diagrams (b) and (c) are described 
by a Brucckncr G-matrix. 



(a) 



(b) 




FIG. 1: Pairing interaction with screening in the RPA ap- 
proximation. The short-dashed line represents the bare in- 
teraction, long-dashed lines the G-matrix, the wiggly line the 
p-h residual interaction resummed to all orders. All vertices 
are to be understood as anti-symmetrised matrix elements. 
The notations are: k = (k, a, r), k — (—k, ~ct, t). 

Within the Brueckner theory of nuclear matter, the 
single particle (sp) energy spectrum of the non super- 
fluid state is derived from the hole-line expansion of the 
mass operator. The BHF approximation, extended to 
include depletion of the Fermi surface due to the strong 
ground state correlations, gives an important quenching 
of the pairing gap and hence it can not be neglected. 
The pp correlations on the pairing interaction are em- 
bodied in the gap equation itself Consistently, the 
interaction vertices in the screening term should be de- 
scribed in terms of the G-matrix. In addition, the lat- 



ter must be dressed according to the Babu-Brown the- 
ory of the induced interaction 0, 0] to avoid the low 
density instability problem of nuclear matter, as dis- 
cussed later. Since the exact resummation of the bub- 
ble series (Bcthc-Salpeter equation) with G-matrix is a 
prohibitive task, the ph vertex insertions can be conve- 
niently evaluated in the Landau limit, and eventually 
replaced by the Landau parameters. 

The effects of the selfenergy corrections have intensely 
been studied; in particular the depletion of the Fermi 
surface is expected to hinder the virtual transitions 
around the Fermi surface and thus its effect is to weaken 
the pairing correlations. It is of particular interest to 
understand to what extent the resulting quenching of 
pairing gap is compensated by in medium vertex correc- 
tions that in nuclear matter are strongly attractive. On 
the other hand, in the case of neutron matter selfenergy 
effects enhance quenching due to screening at variance 
with the predictions of recent Monte Carlo many body 
calculations |l8| . 

Going beyond the pure BHF approximation, the main 
dispersive corrections arise from energy dependence of 
the self-energy, as shown in previous papers. The correc- 
tion to the non superfluid propagator G k (u>) is simply 
a renormalization factor of its pole part. This factor, 
named Z-factor, is 
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which measures the discontinuity of the occupation 
probability around the Fermi energy. Correspondingly, 
the abnormal propagator appearing in the gap equation 
(Eq. (1)) 



A fc (w) 



G7\u)G7\-u) + M{uy 



(3) 



is renormalized by a factor Z 2 . For a static interaction 
the gap function is also independent of energy and the 
energy integration can be performed analytically. Since 
the analytical structure of the abnormal propagator is 
not modified, the gap equation takes the same form as 
in the pure BCS case. One easily obtains: 



d 3 k'V, 



Z k Z k /A k ' 



k.k' 



(4) 



where sf is the Fermi energy and e k is the on-shell self- 
energy. The preceding gap equation is equivalent to the 
BCS version except for the Z 2 factor which is modeling 
the effect of the interaction around the Fermi surface. 
Since the value of the Z-factor, though depending on the 
ground state correlations, is always less than unity in the 
vicinity of the Fermi surface, inevitably the energy gap 
will turn out quenched in this respect. Our predictions 
of the gap in nuclear and neutron matter presented in 
this paper rely on the solution of the latter equation. 
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III. SCREENING INTERACTION 

A. One-bubble screening interaction 

In a previous work the calculation of the screening in- 
teraction was simplified by using the Gogny force, which 
in fact reproduces most of the properties of a G- matrix. 



In the present calculation we adopt the G-matrix itself 
and we try to reduce its complexity with reasonable ap- 
proximations. For the sake of application to the pairing 
in the 1 Sq channel we select the two particle state with 
total spin S=0 and isospin T=l. Then the one-bubble 
interaction can be written as 



J 



< ll|Vi|lT > 



2,2' ST 



-f(2S+l) < l2\G p s h T \l'2 



> A < 2'1|G^|21' > A A°(22') 



(5) 



r 



where 1 = (fci,cri,Ti) (1' = (ky, ay, Ty)) and 1 = 
(— fci,eri,Ti) (I' = (— ky, <Ty,Ty)) are the momenta of 
the pair in the entrance (exit) channel. A is the static 
polarization part. The G-matrix is converted into the 



ph sector, as it is required to solve the Bethe-Salpeter 
equation and to sum up the bubble series V% . The stan- 
dard recoupling procedure from pp sector to ph sector 
yields 



J 



G st = J2(2S C + 1)(2T C + 1)(-1) S °+ T * 
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where the brackets are the 6j symbols. The sum runs 
over the spin S c and isospin T c of the pp channels in- 
cluded in the calculation. Since the G-matrix incorpo- 
rates short range pp correlations, its momentum range 
is shrunk remarkably in comparison with the bare in- 
teraction, as shown in Fig. [21 At variance with the bare 
interaction, the the G-matrix cannot sustain large mo- 
mentum transfers q '— k — k' , that justifies the approx- 
imation to average it around the Fermi surface, in the 
limit q=0. As a consequence the q dependence is only 
located in the integral of the polarization part, giving 
the Lindhard function 



Z ✓ k.k—q 



N(Q) 1 



-1 



1 o 2 

-(l- q T )ln 
q 4 



1-9/2 



l + q/2 



,(7) 



where g is the degeneracy parameter. 

The two external vertices mixing pp and ph lines 
shown in Fig. 1 (b) in principle induce unlimited ex- 
citations in momentum space. But, since high mo- 
menta transitions are incorporated into G-matrix used 
as vertex interaction, the main contribution of the two- 
bubble diagram is concentrated in a domain as short as 
2 fm . Elsewhere the remaining contribution can ne- 
glected with respect to the bare interaction as shown in 
Fig. 2. 

The problems with the calculation of the ph multi 
bubble contribution (diagram (c) in Fig. 1) are the 
following. First, the bubble series with G-matrix in- 
sertions have to be previously summed up. But, since 
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FIG. 2: Shrinking of G-matrix in momentum space. 



the interaction vertices in the ph channel involve par- 
ticle excitations around the Fermi surface, they can be 
approximated by the Landau parameters. Second, even 
replacing the bare interaction vertices by G-matrices , 
there appears the long low density singularity in the 
RPA in nuclear matter (F = —1). This problem, dis- 
cussed in Ref. (see also references therein) is reme- 
died by dressing the vertex insertions according to the 
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Babu-Brown induced interaction theory [3. 



B. Landau parameters from the BHF 
approximation 

The microscopic basis of the ph effective interaction 
can be set in terms of the energy functional of symmetric 
nuclear matter 



N{0)f aT ^ T ,{k,ki) = 



= F + F'(t ■ t') + G(a ■ a') + G'(a ■ ct')(t • r'), 



(8) 



where the density of states ^V(O) is introduced to make 
the Landau parameters F, F' , G and G" dimensionless. 
In BHF approximation the energy functional is given 
by 



E = 



^ 2m 2 

k ki,k2 



^2{k 1 ,k 2 \G{co)\k 1 ,k 2 ) A , (9) 



where the subscript A means that the matrix element of 
the G-matrix is antisymmetrized. The index k stands 
for k, a and r, momentum, spin and isospin, respec- 
tively. The G-matrix is understood to be calculated on 
the energy shell: lu = +£fc 2 - The single particle ener- 
gies are determined iteratively along with the G-matrix 
within the Brueckner selfconsistent scheme. One can 
determine the Landau parameters from the microscopic 
Brueckner theory in the BHF approximation, perform- 
ing the double variational derivative, Eq. (8), of the 
energy per particle, Eq. (9). So doing, a number of 
contributions are generated that can be calculated one 
by one |19| in some approximation due to the complex 
structure of G-matrix. A simple and powerful way to 
calculate the Landau parameters is to suitably fit the 
BHF energy and the corresponding sp spectrum with a 
functional of the occupation numbers and then to per- 
form the double derivative. A Skyrme-like functional 
has proved to reproduce accurately the equation of state 
(EoS) of symmetric as well as spin and isospin asym- 
metric nuclear matter pfjj . Therefore we determine the 
Landau parameters in that way. A limitation of this 
procedure is that only a few partial wave components 
can be calculated, but for the purpose of the present 
investigation we only need the zero order Landau pa- 
rameters. The latter are plotted in Fig.|3|as a function 
of the Fermi momentum. As expected Fq exhibits the 
well known instability below the saturation point, which 
makes the RPA series difficult to handle. As in previous 
papers 0, 0] this drawback can be overcome by the in- 
duced interaction theory of Babu and Brown [la. Leav- 
ing aside a description of this theory (see Refs. [itH^ K 
we schematically write down the equation defining the 
ph induced interaction as follows 



V ph = V d + V RPA (V ph ). 



(10) 



The first term (direct term) is the BHF ph residual in- 
teraction, which, in the first order, is represented by the 



G-matrix. The second one (induced term) is the RPA 
bubble summation, in which the vertex insertions are 
given by V p h itself instead of the direct term. The solu- 
tion of the latter equation is quite simple if we replace 
the true Vd projection in the ph channel (ST) with the 
corresponding Landau parameter, since the way we ex- 
tract the Landau parameters the direct term contains 
not only the effect of the G-matrix but also the rear- 
rangement diagrams |22l | . The numerical results are de- 
picted in Fig. |31 The salient feature of the induced in- 
teraction is that the renormalization of Fq prevents any 
singular behavior to occur below the saturation den- 
sity. Otherwise the values do not differ from the Gogny 
interaction except for the high density behavior where 
the effect of the three body force makes Fq much more 
repulsive than the Gogny force 0. 

Therefore we dressed the residual interaction first 
with the short range correlations (G-matrix instead of 
bare interaction) and then by the renormalized long 
range correlations V p h replacing the G-matrix in the 
RPA series. Since the calculation of the induced interac- 
tion with the G-matrix is a quite complex job, we have 
simplified the problem replacing the G-matrix with the 
Landau parameters. The way we determine the Landau 
parameters, the approximation turns out to be better 
than starting from the G-matrix itself. 
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FIG. 3: Landau parameters of pure neutron matter and nu- 
clear matter. 



C. Bubble series 

The vertex insertions dressing the bubbles must be 
treated on different footing than the external ones (di- 
agram (c) of Fig. 1). To sum up the RPA bubble series 
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with the G-matrix requires to solve the Bethe-Salpeter 
equation, what is a prohibitive task. Using the Lan- 
dau parameters corresponds to the Landau limit (zero 
momentum-energy transfer around the Fermi surface), 
which is a quite reasonable approximation. In this case 
the RPA summation of the ph interaction turns out to 
be algebraic and, expressed in term of the dressed bub- 
ble, it is written as 



A(g) S T = 



A°(g) 



l + A°(q)£ 



ST 



(11) 



where Cst are the Landau parameters, whose compo- 
nents are commonly denoted by: Cqo = F, Cqi = F 1 , 
£10 = G, £11 = G'. In this expression we clearly see 
how the induced interaction prevents any divergence to 
occur since |A£| < \C\ < 1. Replacing in Eq. J5J the 
bare bubble A with the dressed bubble A we get the 
full screening interaction used in the calculation. 



IV. RESULTS 

The G-matrix is generated from a selfconsistent BHF 
calculation with the continuous choice |7j. The Ar- 
gonne V18 two body force is adopted as the input 
bare interaction plus a microscopic model for the three 
body force based on meson exchange with intermedi- 
ate excitation of nucleon resonances (Delta, Roper, and 
nucleon-antinucleon) ■ The calculation also provides 
the self-energy from which we extract the sp spectrum 
and the Z-factors. Based on the same framework is also 
the Skyrme-like fit of the BHF energy functional used 
to calculate the Landau parameters p0| . 



A. Screening interaction 

In this paper we only focus on the 1 Sq pairing interac- 
tion in the two extreme situations of pure neutron mat- 



ter and symmetric nuclear matter. We keep for a fur- 
ther investigation the consideration of asymmetric nu- 
clear matter with the purpose of studying the transition 
from the screening regime in pure neutron matter to the 
antiscreening regime in symmetric nuclear matter. 

Let us start with neutron matter. In this case the 
screening interaction can be decomposed in two terms: 
S=0 density fluctuation and S=l spin density fluctua- 
tion. The two modes have opposite effect: the former 
one is attractive, the latter is repulsive. From Eq. (JSJ), 
we can write 

Vx = l -K{ q ),GfGt - \K{q)iG? G? . (12) 

The factor 3 is due to the multiplicity of the spin mode. 
For the following discussion we should notice that A is 
negative. In turn, each G-matrix can be expressed as a 
superposition of G-matrices, projected onto two particle 
states (see Eq. (6)), 



Gq 11 — 2 (~Gq - 3Gi), 



1 



(13) 
(14) 



Since Go is attractive and G\ is repulsive, assuming 
their magnitude to be comparable, we get Go ~ G\ |24j 
and the multiplicity plays the main role in establishing 
the dominance of the spin density mode over the den- 
sity mode. In the latter calculation Go and G\ are only 
roughly comparable, as it can bee seen in Fig. 4, never- 
theless the conclusion is still valid. A quenching of 1 So 
pairing in neutron matter is to be expected, a result 
established long ago p5l | and confirmed by many calcu- 
lations in various approximations (see and references 
therein). 

In nuclear matter the situation could be quite differ- 
ent since the isospin fluctuations also come into play. 
The screening interaction now is split according to 
Eq. JSJ) as follows 



Vi = \{k(qUGllG§ + A( ? ) 01 GM) - \{k{q) w GtGt + A( ? ) U G#G#). 

I 



The various contributions are plotted in Fig. 4 in terms 
of pp states, the individual ph contributions are ex- 
pressed as 

G P oo = j(G o + 3Gio + 3Goi+9Gh), (15) 
Giq = -(— Goo — 3Gio + G01 + 3Gn), (16) 
Gqi = —(— Goo — 3Gio + G01 + 3Gn), (17) 
G11 = ^-(Gqo — G10 — Gqi + G11). (18) 



In nuclear matter the pp G-matrix elements are domi- 
nated by the deuteron channel ( 3 SDi coupled pp chan- 
nel), which is very attractive and therefore it reinforces 
the density mode and weakens the spin mode. In other 
words, the main isospin effect is to reverse the role of the 
medium, i.e. antiscreening instead of screening. In pre- 
vious papers this effect has been discussed in terms of 
proton-proton ph screening against neutron-neutron ph 
screening in the neutron-neutron 1 So channel [2(j . The 
latter gives repulsion the former attraction. At vari- 
ance with Ref. 24] the proton-proton ph screening is 
stronger than neutron-neutron ph screening. This effect 
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is to be traced back to stronger in medium renormaliza- 
tion of the force in the T — channel than in the T = 
one. Antiscreening is the overall effect. 

In Fig. 5 we plot the full pairing interaction in the 
three approximations used in the calculation of the en- 
ergy gap. In nuclear matter, as we discussed before, the 
screening effects in fact reinforce the attractive strength 
of the bare interaction. The main effect appears al- 
ready at the one bubble level. The deviation from the 
bare interaction increases at lower density. At kp = 0.6 
fm^ 1 the enhancement is from -13 MeV ■ fm 3 to -27 
MeV ■ fm 3 . This is a huge variation which could entail 
a large increase of the gap because it is exponentially 
depending on the interaction. But at such a density 
the pair correlations are rather weak and thus we do 
not expect any large increase of the gap. In the den- 
sity domain of the maximum gap the enhancement is 
much smaller and again we do not expect any dramatic 
change in the gap magnitude as an effect of the anti- 
screening. In Ref. 0, an improper coupling of the ph 
states in the mixed representation prevented the cance- 
lation among different ph excitations to occur with the 
effect of producing a more pronounced antiscreening. 

In neutron matter the situation is the other way 
round. The screening is repulsive, and small in the full 
RPA calculation, but still enough to produce a sizeable 
quenching of the pairing gap. These predictions confirm 
at least at qualitative level the corresponding results ob- 
tained with the Gogny force 0. 
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FIG. 4: Individual components of the ph residual interac- 
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B. Pairing gap 

The present calculation is focussed on the 1 5'o 
neutron-neutron (or proton-proton) pairing. One can 
distinguish the bare interaction which is responsible for 
the pairing between the two particles in the 1 Sq state, 
from the screening interaction induced by the surround- 
ing particles. Therefore the interaction, projected onto 




K F (fm" 1 ) 

FIG. 5: Pairing interaction. 

the 1 Sq channel, can be cast as follows 

<k\V\k'>= / ^[VoCMO+Vidfc-fc'l)]. (19) 

As bare force we use Argonne V18, the same as for 
calculating the G-matrix and the selfenergy. The BCS 
energy gap in the 1 So channel is practically independent 
of the adopted bare force, since in fact all realistic inter- 
actions reproduce the phase shifts of free NN scattering. 

We solved the gap equation in the form of Eq. (0J. 
In order to disentangle the screening effects from the 
selfenergy corrections, we first assume Z—l and free sp 
spectrum. The results are plotted in Fig. (upper left 
panel). In neutron matter the screening effect is small 
and just reduces the gap by 10% in the peak region. At 
variance with previous calculations existing in the litera- 
ture the full RPA screening is much less effective than 
the one bubble approximation because of the stronger 
renormalization of the spin fluctuations vs the density 
fluctuations in the induced interaction. However this 
finding confirms the preceding predictions with Gogny 
force (see Fig.8 of Ref. B). 

In nuclear matter, due to the antiscreening effect we 
discussed earlier, the magnitude of the gap variation is 
the other way around and much more sizeable: the gap 
rises up from 3 M eV to 5 MeV for Fermi momentum 
kp = 0.8 fm -1 . This is displayed in Fig. 6 (lower left 
panel). 

There are two kinds of selfenergy effects: dispersive 
effect and Fermi surface depletion. Both are calculated 
taking into account the selfenergy corrections at the sec- 
ond order of G-matrix (rearrangement terms). The first 
one is a correction to the sp spectrum in the energy de- 
nominator. Usually it entails a reduction of the pairing 
gap since the effective mass, beyond BHF approxima- 
tion, is less than the unity (the effective mass is the 
combination of the e-mass and the k-mass Q). But at 
very low density the effective mass is larger than unity 
|27| and it reduces the quenching rate of the gap due 
to the interaction. This effect can be seen in the low 
density side of the neutron gap with Ti to tai (upper right 
panel). Additional strong reduction is due to the de- 
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plction of the Fermi surface which hinders transitions 
around the Fermi surface. The maximum gap in a com- 
plete calculation is 1.5 — 2 MeV at hp w 0.8 fwT 1 . 

In nuclear matter the selfenergy effects are much 
stronger already at moderately low density, as it has 
to be expected, and the peak value shifts down to very 
low density k,F ~ 0.5 — 0.6 /to -1 . The Z-factor plays 
the major role: it quenches from 0.84 in neutron matter 
to 0.68 in nuclear matter at Uf = 0.8 /to -1 . But the 
magnitude is about 0.5 MeV less than the value with 
only bare interaction. Therefore we can conclude that a 
strong cancelation occurs as soon as vertex corrections 
and self energy effects are simultaneously included in the 
gap equation. But this happens only in nuclear matter 
as an effect of antiscreening. We will come back to this 
point below. 
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V. DISCUSSION AND CONCLUSIONS 

In this paper an exhaustive treatment of the 1 So pair- 
ing in nuclear and neutron matter has been reported. 
The medium polarization effects on the interaction and 
the selfenergy corrections to the mean field, both devel- 
oped in the framework of the Brueckner theory, have 
been included in the solution of the gap equation. 

Within the pure mean field approximation [2^| the 
1 Sq gap is not affected by the medium, either nuclear 
or neutron matter. So far the medium effects have not 
been considered in the case of nuclear matter except in 
Ref. The vertex corrections due to neutron matter 
all give a reduction of the pairing, the magnitude de- 
pending on the adopted approximation [3j- The expla- 



nation relies on the competition of the attractive den- 
sity excitations against the repulsive spin density excita- 
tions. The present calculation, based on G- matrix, also 
predicts a large quenching in agreement with almost all 
previous predictions, but only at the one bubble level. 
In the most complete calculation (full RPA) the quench- 
ing is largely reduced in apparent agreement with a re- 
cent Monte Carlo calculation But the inclusion of 
selfenergy effects definitely results in a large suppression 
as expected from basic properties of a strongly corre- 
lated many body system (see Introduction). 

In the case of nuclear matter the most remarkable re- 
sult is the antiscreening effect of the medium polariza- 
tion. In fact in nuclear matter isospin modes arise that 
reverts the competition between the attractive density 
modes and the repulsive spin-density modes due to the 
presence of isospin modes. The argument addressed in 
Ref. [24| that the p-n (T=0) interaction is small com- 
pared to the n-n (T=l) is based on the vacuum scatter- 
ing T-matrix and does not consider the strong medium 
renormalization of G-matrix, which inverts the strength 
of the two channels. However the enhancement of the 
gap to almost 5 MeV is almost completely suppressed 
by the strong correlation effects on the selfenergy. But, 
even a small variation of the force strength implies a 
large variation of the gap. These effects also push to 
lower density the peak value of the gap. 

Calculations of the pairing gap in a nuclear environ- 
ment have been reported in a series of papers for the 
case of nuclei Their main finding is that the induced 
interaction arising from the surface vibrations is respon- 
sible for large part of the experimental gap. This result 
can be considered as the counterpart for finite nuclei of 
the antiscreening effect due to the medium polarization 
in nuclear matter. But the selfenergy effects completely 
compensate the gap enhancement and, in the end, the 
full medium effect do not change significantly the gap 
with bare interaction. This result turns out to be not a 
big surprise, since some calculations show that the gap 
with Gogny interaction is consistent with the observed 
gaps in nuclei [29j| . 

At this point two aspect are worth to be developed 
further. The first one is the study of pairing in the 
transition from symmetric nuclear matter to neutron 
matter; the second one is the investigation of pp chan- 
nels so far neglected, since there the bare interaction is 
repulsive. The pairing could exist as an induced effect 
of the environment. 
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